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Distinguishable and indistinguishable particles 



In the previous parts of this course, we derived the Boltzmann 
distribution, which described how the number of distinguishable 
particles in different energy states varied with the energy of 
those states, at different temperatures: 



However, in systems consisting of collections of identical 
fermions or identical bosons, the wave function of the system 
has to be either antisymmetric (for fermions) or symmetric (for 
bosons) under interchange of any two particles. With the 
allowed wave functions, it is no longer possible to identify a 
particular particle with a particular energy state. Instead, all 
the particles are "shared" between the occupied states. The 
particles are said to be indistinguishable. 
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Indistinguishable fermions 



In the case of indistinguishable fermions, the wave function for 
the overall system must be antisymmetric under the 
interchange of any two particles. One consequence of this is 
the Pauli exclusion principle: any given state can be occupied 
by at most one particle (though we can't say which particle!) 

For example, for a system of two fermions, a possible wave 
function might be: 

V>Oi, x 2 ) = A= [^a(xi)^ b (x 2 ) - ^a(x 2 )^b(xi)] . (2) 



Here, x\ and x 2 are the coordinates of the two particles, and A 
and B are the two occupied states. If we try to put the two 
particles into the same state, then the wave function vanishes. 



Statistical Physics 2 Part 5: The Bose-Einstein Distribution 



Indistinguishable fermions 



Finding the distribution with the maximum number of 
microstates for a system of identical fermions leads to the 
Fermi-Dirac distribution: 



9i 

rii = • (3) 

e kT + 1 



Here, is the number of particles with energy e^; is the 
degeneracy of the energy level (i.e. the number of distinct 
states with energy ej; k is Boltzmann's constant; and T is the 
thermodynamic temperature. The function /x(T) is known as 
the chemical potential, and is determined from the condition: 

T, n i = N > ( 4 ) 

i 

where N is the total number of fermions in the system. 
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Indistinguishable fermions 



Since e x > 0 (for real values of the variable x), we find that for 
all values of e{ and T: 

Tin 1 , . 

i = ^ — ^- ( 5 ) 

9% e ~w~ + 1 



This tells us that the number of fermions with energy ^ must 
be less than or equal to the number of states gi with that 
energy. This is another statement of the Pauli exclusion 
principle. 
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Indistinguishable bosons 



The wave function for a system of identical bosons must be 
symmetric with respect to the interchange of any two of the 
bosons. For example, with two bosons, a possible wave 
function might be: 

tp(x 1 ,x 2 ) = [VmOi)^^) + ^a(x 2 )^b(xi)] . (6) 

This time, if we put two particles into the same state, the wave 
function does not vanish: there is no limit on the number of 
particles we can put into any given state. 

The Pauli exclusion principle does not apply to bosons: and 
this implies that the Fermi-Dirac distribution does not apply 
either, since the Pauli exclusion principle follows from the 
Fermi-Dirac distribution. 
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The Bose-Einstein distribution 

In this part of the lecture course, we shall derive the 
Bose-Einstein distribution for a system of indistinguishable 
bosons, and apply the distribution to three cases: 

• the Bose-Einstein gas; 

• a collection of photons; 

• the vibrations of particles in a solid. 



To find a distribution function to describe systems of bosons, 
we shall follow the same procedure that we did for 
distinguishable particles, and for indistinguishable fermions. 
That is, we write down the number of microstates for the 
possible distributions, subject to constraints determined by the 
macrostate specifications. Then we look for the distribution 
with the largest number of microstates: applying the principle 
of equal a priori probabilities, this distribution is more likely to 
occur than any other. 
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The Bose-Einstein distribution 



Consider an energy level ^ with degeneracy containing 
bosons. The states may be represented by gi — 1 lines, and the 
bosons by circles; distinguishable microstates correspond to 
different orderings of the lines and circles. For example, with 9 
particles in 8 states corresponding to a particular energy, a 
particular microstate might be: 



The number of distinct orderings of lines and circles is: 

(V>i + gi - 1)! 



k = 



riiKgi-iy. 



(7) 
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The Bose-Einstein distribution 



A particular distribution has a specified number of particles ni 
within each of the possible energy levels e^. The total number 
of microstates for a given distribution is therefore: 

U /; V ni!(^-l)! 



Let us assume that each state has a high degeneracy, i.e. 
gi > 1. Then we can make the approximation: 

«({n,})*n ( "'t 9 | )! - ( 9 > 
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The Bose-Einstein distribution 



To find the most probable distribution, we follow the same 
procedure as for the Boltzmann and the Fermi-Dirac 
distributions. That is, we maximise (9): 

Oz + 9iV- 



= n 



, ni\gi\ 



subject to the constraint on the total number of particles: 

T, n i = N > ( 10 ) 

i 

and the constraint on the total energy: 

5><n< = tf (11) 

i 

As usual, rather than maximise t directly, we maximise Int. If 
we assume that both gi and are large enough for Stirling's 
approximation to hold for ln^! and Inn^!, we find that ln£ is 
given by: 

In t « K n i + 9i) ^ (ni + 9i)-9i\ngi- In n 4 ] . (12) 

i 
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The Bose-Einstein distribution 



The change in ln£ resulting from changes dni in each of the 
populations ni is then: 

d In t « ^2 [' n ( n z + 9i) — In drij\ . (13) 
From the constraints (10) and (11), we find: 

Y^dni = 0, J2 £ i dn i = °- ( 14 ) 

% i 

Combining (13) and (14) with Lagrange multipliers a and (3, 
we have: 



din t « 



I n + a + pe 



dni. (15) 



For appropriate values of a and (3, equation (15) is true for all 
dni, hence: 

ln (^^) +CK + / ^ = °- (16) 
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The Bose-Einstein distribution 



We then find that the most probable distribution can be 
written: 

ua — § . (17) 

1 p -a-0Si _ i 



Equation (17) is the Bose-Einstein disbribution. It gives the 
population of an energy level that has energy ^ and degeneracy 
flfj. The constants a and (3 are determined from the constraints 
(10) and (11) on the total number of particles and the total 
energy. (3 can, as usual, be related to the thermodynamic 
temperature, so that the Bose-Einstein distribution takes the 
form: 

ni = • (18) 

BekT - 1 
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Three distributions 



We have now completed our collection of three distributions. 
First, a collection of distinguishable particles follows the 
Boltzmann distribution: 

N 

nj = —e kT. (19) 
Note that j indexes particle states, not energy levels. 

Then, for a collection of indistinguishable fermions, we found 
the Fermi-Dirac distribution: 

"> = -T^^- (20) 

e kT + 1 

In this case, i indexes energy levels, not individual states. 

Finally, we found the Bose-Einstein distribution for a collection 
of indistinguishable bosons (18): 

9i 



Hi = 



BekT - 1 
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Fermions and bosons 



The formulas for the Fermi-Dirac distribution (20): 

— 9i 

n i — ej-fi 

e kT + 1 

and the Bose-Einstein distribution (18): 

9i 



Tin = 



BekT - 1 

appear very similar. However, the difference in sign in the 
denominator is critical. Since e x > 0 for all x, for the 
Fermi-Dirac distribution, we always have < gi\ in other 
words, fermions always obey the Pauli exclusion principle. 

But for bosons, it is possible to have > g it if ^ <c kT, and 
B«l. In other words, it is possible for many identical bosons 
to exist in the same state. 
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Bose-Einstein distribution: low-density limit 

When we derived the Fermi-Dirac distribution, we showed that 
for systems with low density, the distribution could be closely 
approximated by the Boltzmann distribution. Using the same 
reasoning, we can show that this is also the case for the 
Bose-Einstein distribution. 

We write the Bose-Einstein distribution (18) in the form: 

V± — 1 

9i BekT - 1 

If rii <c gi, then rii/gi <c 1, which implies that: 

BekT - 1 > 1, 

so we can approximate: 

BekT - 1 « BekT. 

We then have: 

9i -^L 
ns K,—e kT. 

1 B 
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(21) 

(22) 

(23) 
(24) 



The Bose-Einstein gas 



As an example of the Bose-Einstein distribution, let us consider 
a boson gas. This consists of a large number of identical 
bosons in a box with rigid walls and fixed volume. The bosons 
are free to move within the box, but cannot move beyond its 
walls. Bosons exchange energy in occasional collisions with 
each other. For simplicity, we assume that the bosons have 
spin zero, so there is no degeneracy associated with different 
spin states. 

Let us first consider the behaviour of the constant B that 
appears in the expression for the Bose-Einstein distribution. 
This quantity corresponds to the partition function in the 
Boltzmann distribution, or the chemical potential in the 
Fermi-Dirac distribution. B is determined by the constraint: 

T, n i = N > (25) 

i 

where N is the total number of particles. Let us find how B 
depends on temperature. 
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The Bose-Einstein gas 



From our previous analysis of the Maxwell-Boltzmann gas and 
the Fermi gas, the density of states for spin-0 particles moving 
freely in a box of volume V is: 

V /2m\§ n 
9(£) = 4^(^J ^ (26) 
Let us assume that we can replace a summation over discrete 
energy levels by an integration over a continuum of energy 
levels. Then the constraint (25) becomes: 

J2 n i~ n(e)de = N. (27) 

Using the Bose-Einstein distribution (18) and the density of 
states (26), this becomes: 



3 ,- 

V [2m\2 r<x> y/e 



4tt 2 V h 2 J JO Be^T - 1 



-de = N. (28) 
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The Bose-Einstein gas 



To simplify things a little, we define a new variable, y, such 
that: 



y = w (29) 



In terms of y, the equation for B becomes: 

3 



47T 

We define the function F(B) 



V /2mAT\ 2 roo ^/ydy_ _ 

-7T 2 V h 2 J Jo BeV - 1 V J 



F{B) = ^r^L. (31) 

Vtt Jo BeV — 1 



In terms of F(B), the equation for B becomes: 



3 

(mkT\i , x . « 
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The Bose-Einstein gas 



To understand the function F(B), let us draw a graph of it: 
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The Bose-Einstein gas 



For large B, we see that F(B) » 1/B. This is because, for 
B > 1, we can write: 

2 foo ^/ydy 2 foo ^/ydy _ 1 



V^r" Jo BeV - 1 Jo BeV B 



(33) 



Thus, from equation (32): 



,'mkT\2 , N 



= iV, 



we find that, for large B (or, equivalently, for high 
temperature): 



B 



V /mkT\2 



N \2irh 2 



(34) 
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The Bose-Einstein gas 



However, we find that the integral cannot be evaluated below 
B = 1, because there is a singularity at this point. Physically, 
this is significant. The Bose-Einstein distribution (18) is: 

9i 



m = 



BekT 



If B < 1, there is the possibility of an energy level ^ such that 
rii < 0. Clearly, we cannot have a negative number of particles 
in any energy level, so we expect B > 1. 
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The Bose-Einstein gas 



We see that F(B) takes a maximum value, at B = 1. The 
maximum value of F{B) is given by: 

W V^Jo eV - 1 S V27 V J 



Here, ((x) is the Riemann zeta function, defined by: 

CO*) = E (36) 
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The Bose-Einstein gas 



The fact that F(B) has a maximum value of « 2.612 is 
somewhat disturbing. It is disturbing because, since N is 
constant, we expect from equation (32): 

that as the temperature goes down, F{B) has to go up. 



But if F(B) has a finite maximum, then below some 
temperature T B , the number of particles we find by integrating 
the Bose-Einstein distribution over all energies is less than N, 
the number of particles in the system: 

V f2mkT\i roo x /ydy , T ^ 
ifT<T B , then — = ( — =— ] / v < N for a B. 
B: 4tt 2 \ h 2 Jo BeV - 1 
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The Bose-Einstein gas 

The situation becomes clear if we plot the population density 
n{e) as a function of energy e, for several different 
temperatures. 

For T > T B , we see a smooth distribution curve that, for high 
temperatures, approaches the distribution expected for a 
Maxwell- Boltzmann gas. The area under the curves for T > T B 
are independent of the temperature T. 

For T < T B , strictly speaking we cannot evaluate the population 
density, since we don't have a solution for the parameter B. If 
we do the best we can, and set B = 1 for T < T B , we find a 
distribution that is sharply peaked at low energy: but the peak 
is not high enough for the area under the curve to give the 
correct result for the number of particles in the system. 
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The Bose-Einstein gas population distribution 
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Reminder: the Fermi gas population distribution 



Compare the population distribution for the Bose-Einstein gas 
with that for a Fermi gas. The population distribution for a 
Fermi gas is well-behaved even in the limit of low temperature, 
T -> 0. 
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The Bose-Einstein condensation 



The temperature T B at which particles start to "disappear" is 
given by: 



that is: 



T B = m^j- OS) 



The natural question is, where are the "missing" particles for 

T < T B ? 

The answer is: in the ground state! 
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The Bose-Einstein condensation 

The problem, which we now need to rectify, is that when we 
replaced the summation over discrete energy levels by an 
integral over a continuum of energies in equation (27): 

l^ni « n(e)de = N. 
we omitted the particles in the ground state. 

Since the lower limit of the integral is zero, we should assume 
(for consistency) that the ground state has zero energy. But 
the distribution function g{e) oc y/e; so at zero energy, the 
density of states is zero, and the population (in the small 
energy range 0 to de) is also zero. 

In other words, the integral in equation (27) only counts 
particles in the excited states, and ignores those in the ground 
state. 
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The Bose-Einstein condensation 



We made the same mistake when we considered the Fermi gas: 
but in that case, there could only be two particles in the 
ground state (assuming degeneracy two, for spin-^ particles). 
In other words, we omitted two particles out of a large number 
N, and this did not affect our conclusions. But for a 
Bose-Einstein gas, there is no limit to the number of particles 
that can fall into the ground state. So if we omit the ground 
state from the sum, we risk making a large error. 

Our reasoning implies that at temperatures T > T B , all particles 
are in excited states, and the ground state is essentially 
unoccupied. But at temperatures T < T B , particles gradually 
fall into the ground state, which eventually, in the limit T^O, 
contains all the particles in the system. 
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The Bose-Einstein condensation 



The particles in the ground state for T <T B constitute a 
Bose-Einstein condensation. The temperature T B is known as 
the Bose temperature. 

Below a critical temperature T B in a system of indistinguishable 
bosons, the population of the ground state in a series of 
quantised states becomes very large, and discontinuous with 
the populations of the excited states. The extra population of 
the ground state is the Bose-Einstein condensation. 

The existence of the Bose-Einstein condensation leads to 
dramatic effects in systems of bosons at temperatures below 
the Bose temperature. These are discussed further in the part 
of this course covering low-temperature physics. 
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The Bose-Einstein gas 



The population density in the excited states (i.e. excluding the 
ground state) in a Bose-Einstein gas is given by: 



n e x(e) — 



V f1m\ 2 ^/e 



4tt 2 V h 2 J BeXT - 1 



(39) 



The function B{T) plays the role for the Bose-Einstein 
distribution that the chemical potential did for the Fermi-Dirac 
distribution. For a Bose-Einstein gas, B(T) is given by: 



B(T) = 



r i 

F' 1 



2) V T J 



for T < T B , 
for T > T B . 



(40) 
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The Bose-Einstein gas 



The function F 1 is the inverse of F, given by (31): 

_ 2 roo ^/ydy 



2. r ( 
F{B) = — 

\/7T JO 



Jo BeV - 1 



The Bose temperature T B is given by (38): 

2tt// 2 / A r 



2 
3 



T B = 



mk 



,c(lW ' 



Note that, since F(B) — > 1/B as B — > oo, we have: 

1 



B 



x 3 

T \ 2 



for T > T B . 



(41) 



Statistical Physics 



31 Part 5: The Bose-Einstein Distribution 



The Bose-Einstein gas 
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The Bose-Einstein gas: excited and ground states population 
The total population N ex in the excited states is, for T < T B : 



N ex = 



V (2m\2 f°° y/ede 



47r 2 V H 2 J JO elk - I 



(42) 



_ V /2mkT\2 roc ^fydy 

Jo 



eV-1 



(43) 



(44) 



= N 



v 3 

T \ 2 



(45) 



This leaves the remaining population TVq in the ground state: 



N 0 = N 



3 

T \ 2 



T 



B, 



for T < T B . 



(46) 
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The Bose-Einstein gas: ground state population 
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The Bose-Einstein gas: total energy 



If we assume that the ground state is a state of zero energy, 
then particles in the Bose-Einstein condensation make no 
contribution to the total energy of the Bose-Einstein gas. The 
total energy U can then be written: 



roc 

U = en ex (e) de — 
Jo 



3 3 

A^VW) Jo BeW-1 
where we have used equation (39) for the density of particles in 
excited states n ex - Changing the variable of integration, the 
total energy can be written: 



3 3 

U = — ^ [ — ) (/cT)2 ' J 



4tt 2 V h 2 



0 Bey -1 



(48) 



Finally, using equation (38) for T B , we can express the total 
energy as: 



U = 



NkT 



3 3 

T \ 2 r°° y 2 dy 



0 BeV 



(49) 
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The Bose-Einstein gas: total energy 



At low temperature (T < T B ), we can take B = 1, and the total 
energy becomes: 



U = r*s N kT 



x 3 3 



_ 3C(§) 



7V/cT 







3 












3 




2 


s) 





r°° y 2 dy 
Jo eV-r 



2 c(§) 

O.UONkT I — ] . (50) 



We see that at low temperature (T < Tg): 



5 



UocT2. (51) 



Statistical Physics 36 Part 5: The Bose-Einstein Distribution 



The Bose-Einstein gas: total energy 
At high temperatures, B > 1, we can make the approximation 

B e y - i ^ B e y-. 



U « T^vNkT 



x 3 3 
T \ 2 1 fooyldy 



T B J B Jo eV 



3 / T \ 2 1 



~ 2C(|) 
Using (41) for B, we find 



JVfcTl-l -. (52) 



U « |iVfcT. (53) 



Note that this is the same result as for the Maxwell- Boltzmann 
gas: as we might expect, since the Boltzmann distribution 
approximates the Bose-Einstein distribution at high 
temperature. 
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The Bose-Einstein gas: total energy 
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The Bose-Einstein gas: heat capacity 



The heat capacity is given, as usual, by the derivative of the 
total energy with respect to temperature: 




(54) 
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A "gas" of photons 



We turn now to the next system we shall analyse using 
Bose-Einstein statistics: a "gas" of photons. Photons have 
spin-1: they are bosons, and their statistics should therefore be 
described by the Bose-Einstein distribution. 

However, we need to make a significant modification to our 
previous work. Our derivation of the Bose-Einstein distribution 
included the constraint (10), that the number of particles in 
the system was fixed: 



This is appropriate for a system of charged, stable (i.e. 
non-decaying) bosons, such as helium-4. However, the number 
of photons in an enclosed space is not conserved: the total 
number of photons can change over time, as photons are 
absorbed by and emitted from the walls of the enclosure. 
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Our previous analysis included the constraint that the total 
number of particles was conserved; in that case, the analysis 
led to the Bose-Einstein distribution (17): 



The constraint that the total number of particles is constant 
can be lifted by setting a = 0. We then find, for a system of 
photons, that the number of photons n{e)de in a (small) 
energy range e to e + de is given by: 



T, n i = N - 



A "gas" of photons 



9% 




- 1 



n(e) de — 



g{e)de 



(55) 



e 



e kT - 1 
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Photon energy distribution 



To complete the expression for the photon energy distribution, 
we need an expression for the density of states g{e). Since 
photons in an enclosed space are subject to the same boundary 
conditions that we assumed in the case of the 
Maxwell- Boltzmann gas, we can simply use the result (from 
Part 3 of this lecture course): 



V 

g (k) dk = 2x 4tt — ^ k 2 dk, (56) 



where V is the volume of the container, and the initial factor 2 
comes from the fact that there are two polarisations for the 
photons. 
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Photon energy distribution 



Now we use the (quantum) relationship between photon energy 
and wave number: 

e = hv = — = he — = hck. (57) 
A 2tt 



Since the energy is simply proportional to the wave number, we 
find: 

V £ 2 

g(e) d£ = 8tt de. (58) 

(27T )3 {hc)^ 



Thus, the number of photons n{e)de in the energy range e to 
£ + d£ becomes: 

, N , V £ 2 d£ 8ttV £ 2 d£ 

n(£)d£ = 8tt — = — -3 — . (59) 

(2^)^ (he) 6 e w - 1 (he) 6 e w - 1 
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Photon energy spectrum 



Using the previous result (59), the energy u(s)ds in the energy 
range e to e + de is given by: 

u(,eja£ = n(£j£ae = - — -~ — g . (60) 

(Jic)-* e feT - 1 



Conventionally, we write the energy density u as a function of 
the photon frequency v. Since the energy and frequency are 
related by the quantum formula: 

e — hv, (61) 
we find that the energy density is: 

8ttV (hv) 3 dv 
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Photon energy spectrum: Planck radiation formula 



The energy density in the radiation as a function of frequency 
and temperature is given by (62): 

8i:V (hv) 3 dv 



u(v) dv — 



This is the Planck radiation formula, giving the energy 
spectrum of electromagnetic radiation in thermal equilibrium at 
temperature T. 

The predictions from the formula are in excellent agreement 
with experimental data. 
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Photon energy spectrum: Planck radiation formula 
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Photon energy spectrum 



It is interesting to consider three features of the energy 
spectrum: 

• the position of the peak in the spectrum; 

• the total energy; 

• the average photon energy. 

We shall consider how each of the above quantities varies with 
temperature. 
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Photon energy spectrum: position of the peak 



Let us consider first the peak in the spectrum. This occurs at: 

— = 0. (63) 

av 

After some algebra, we find that this leads to the result that 
the peak occurs at: 

l--x)e x =l, (64) 



3 

where x = hv/kT. The (numerical) solution gives: 

kT 

9^2.82 — , (65) 
h 

i.e. the frequency at which the energy spectrum peaks is 
proportional to the temperature: this is Wien's law. 
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Photon energy spectrum: total energy 

Now let us consider the total energy, U. This is the area under 
the graph of the energy spectrum, i.e.: 



JO h 2 c 3 Jo M 1 V J 



n c ju ekT _ i 

We define a new variable of integration: 

hv 

y = w (67) 
in terms of which, the expression for the total energy becomes: 

u= w cm *rpL. (68 ) 

{hey Jo eV — 1 

The integral has the value 7r 4 /15, hence: 

u = H? {kT)A - (69) 
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Photon energy spectrum: average photon energy 


The average photon energy is found from: 




1 r°° 
e = — / n(£) £ de, 
N Jo W 


(70) 


WIlclc iv lb Lllc LULdl lIUlilUcl Ul pilULUllb. 




rOC 

N = n(e) d£. 


(71) 


Using equation (59) for the number of photons per unit energy 


range: 




8irV £ 2 d£ 

n(£)d£ = 

{hey e wr - 1 




we find: 




-_ ekT-l 
ekT-l 


(72) 
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Photon energy spectrum: average photon energy 


Again defining a new variable of integration y = e/kT, the 




average photon energy becomes: 




coo y 3 i 
__/0 e»-l* fcT 

Jo e y_i^ e 


(73) 


Using the results for the integrals: 




fOC y 3 7T 4 fOO y 2 

L v — T de = Tr' L -^r<fe = 2C(3)«2.40, 
Jo ey — 1 15 Jo eV — 1 


(74) 


where C is the Riemann zeta function, we find: 




£ = — /cT « 2.70 x /cT. 

30C(3) 


(75) 


We see that the average energy of the photons increases 




linearly with temperature. 
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Black body radiation 



We can use the spectrum of radiation in thermal equilibrium to 
find the radiation power from an object, as a function of the 
temperature of that object. 

Consider an object that completely absorbs all radiation falling 
on it, of whatever wavelength. Such an object will appear 
completely black. If the object is in thermal equilibrium with a 
"gas" of photons at temperature T, then the object must emit 
as much radiation as it absorbs, and with the same spectral 
properties. 

If the object emits either more or less radiation than it absorbs, 
then the energy in the radiation at some wavelength(s) will 
either increase or decrease; in other words, the object cannot 
be in equilibrium with the radiation. 
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Black body radiation 



Therefore, if we calculate the rate at which the object absorbs 
radiation, as a function of temperature, then the same 
expression gives the rate at which the object emits radiation. 
Since we need to assume that the object absorbs all the 
radiation that fall on it, we will calculate the temperature 
dependence of "black body" radiation. 

To calculate the rate at which a black body emits radiation, we 
consider radiation in a box of fixed volume V. The radiation is 
in thermal equilibrium with the box, which absorbs all the 
radiation that falls on it. The rate at which the inside walls 
emit radiation is equal to the rate at which radiation falls on 
the walls. 
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Black body radiation 



Consider a small area element A/1 on an inside wall of the box. 
Take the centre of this area element to be the origin of a 
spherical polar coordinate system. A volume element dV of the 
interior of the box is given by: 



dV = r Sin 9 dr d9 dep. 



(76) 



If the total number of photons in the box is N, and the volume 
of the box is V, then the number of photons in the volume 
element dV will be: 



dV No 

N = — r z S\n 6 dr dO deb. 

V V 



(77) 
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Black body radiation 



rsin^ d<f> 
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Black body radiation 



Assuming that the radiation is isotropic, the photons inside the 
volume element dV will be travelling with equal probability in 
all directions. 

Therefore, the probability that a given photon in the volume 
element dV will be incident on the area element A/1 is the 
effective area of the area element (taking into account its 
orientation with respect to a line joining it to the volume 
element), divided by the surface area of a sphere centered on 
the volume element, and passing through the area element: 

costfAA . . 

Ap = - . (78) 
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Black body radiation 



The number of photons passing through the area element A/1 
in time At is found by integrating y/±p over all volume 
elements within a distance r = cAt of the area element: 



f N cos 0 A/1 . . 

n 7 AAAt = / - dV, (79) 

= / dr cos 6 sin Odd / d6, (80) 

V 4tt Jo Jo Jo 

NAA 1 , N 

= xcAt x-x2?r, (81) 

N 

= cAAAt. (82) 

AV V J 
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Black body radiation 



Using N = U/i, we find that the radiation energy falling on the 
inside walls of the box per unit area per unit time is (the "flux" 
of radiation): 

U , V 

Pj — rij£ — (83) 



4V 



From equation (69) for the total energy: 



the flux of radiation is: 



2tt 5 (/cT) 4 , . 

^ = -15 W' (84) 



Statistical Physics 58 Part 5: The Bose-Einstein Distribution 



Black body radiation 



As we argued above, since the box is in thermal equilibrium 
with the radiation, the rate at which the walls emit radiation 
must equal the rate at which radiation falls on the walls (and is 
absorbed). 

We now claim that a perfectly absorbing surface at some 
temperature T will emit radiation at the same rate, whether or 
not the surface forms the inside of a box. 

Hence, the rate of radiation of energy per unit area from a 
perfectly absorbing surface at temperature T is given by (84): 



7 

where: 



p = a T 4 , (85) 



2 ^ 5 kA 

° = -15W (86) 



This is Stefan's law; a is Stefan's constant. 
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Radiation pressure 



Finally for the photon gas, we note that since the total energy 
U depends on the volume V, we expect the radiation to exert a 
non-zero pressure on the walls of its container. From the first 
law of thermodynamics, we find that: 

'dU\ 



p = -{av)s- (87) 

The total energy can be expressed in terms of the number of 
particles in the various energy states: 

U = J2 n j £ j- ( 88 ) 

j 

The entropy S is related (by Boltzmann's equation) to the 
total number of accessible microstates, which in turn is related 
to the number of particles in each state. Hence, the constraint 
that S is constant is equivalent to the constraint that the 
energy state populations nj are constant: 

3 
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Radiation pressure 



Now, for photons, the energy of a state is proportional to the 
wave number, k: 

de , dk . . 

e — hck, so: — = he — . (90) 

dV dV v y 

The boundary conditions impose constraints on the possible 

wave numbers. In a cube of side length L, we have: 

k = ^^nl + nl + nl, (91) 

for positive integers %, n y and n z . This can be written: 

2tt 



from which it follows that: 



k = + nl + n 2 z , (92) 



dk 1 k de 1 e , 

— = , so: — = . (93) 

dV 3V' dV 3V v y 
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Radiation pressure 



Substituting (93) into (89) gives: 



p= -Vn-^ = -— (94) 
3^ J V 3V' 

3 



The radiation pressure depends on the energy density in the 
radiation. 

Using equation (69) for the total energy gives: 

87r5 ^ ro^ 
p = ^Jh^' (95) 

The radiation pressure depends only on the temperature (and 
some fundamental physical constants); and increases as the 
fourth power of the temperature. 
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Pressure of a gas of massive bosons 



Note that equation (94) for the pressure in a "gas" of photons: 

1U 

P=3Y (Photons) 

differs from the expression we would find for a gas of massive 
bosons. In the latter case, the energy of a state is quadratic in 
the wave number k, rather than linear: 

- = I* 2 - (96) 
2m 



We have the same boundary conditions for massive bosons as 
for massless photons. The result is that the pressure in a gas 
of massive bosons depends on the energy density as: 

2U 

p= — (massive bosons). (97) 
3V 
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Final example: heat capacity of a solid 



In Part 2 of this course, we considered a simple model of a 
solid, in which the atoms were represented as weakly-coupled 
quantum harmonic oscillators. The energies of the 
single-particle energy eigenstates were then: 

£j = (^ + j)hu, (98) 

where j is zero or a positive integer, and a; is a constant, 
characteristic of the force between atoms, and the mass of 
each atom. This is the "Einstein model" of a solid. 



With the above energy levels, we were able to perform an 
analysis of the system using statistical mechanics, to derive 
quantities such as the total energy, entropy and heat capacity 
as functions of temperature. For the heat capacity, we found 
qualitative agreement with experiment, in that the heat 
capacity goes towards zero at zero temperature. However, 
quantitatively, the predictions are not very accurate. 
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Heat capacity of a solid 
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Heat capacity of a solid: Debye model 



A more sophisticated model of a solid has been suggested by 
Debye. In Debye's model, we consider not the vibrations of 
individual particles, but "collective" vibrations, comprising 
standing waves in the solid. The key difference when compared 
to the Einstein model, is that the frequency of a wave depends 
on its wavelength: 

• In the Einstein model, the vibrations have fixed frequency, 
and the energy levels are uniformly distributed in energy. 

• In the Debye model, the vibrations cover a range of 
frequencies, and the density of states varies with frequency. 
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Debye model: cut-off frequency 



Another key feature of the Debye model is the cut-off 
frequency. A solid can only support waves with frequencies 
below a set limit. This is equivalent to saying that waves in the 
solid have a minimum wavelength. The minimum wavelength is 
determined by the spacing of atoms in the solid: clearly, there 
cannot be waves with wavelength less than twice the spacing 
between atoms. 

To find the cut-off frequency, let us consider standing waves in 
the solid. We shall assume the boundary condition that atoms 
at the surface of the solid are in fixed positions (i.e. atoms at 
the surface cannot vibrate). To make things simple, let 
consider initially a one-dimensional solid consisting of just a 
small number of atoms... 
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Debye model: vibrational modes 



Let the atoms in the one-dimensional solid be numbered 
m = 0, 1, 2, N. The total number of atoms is N + 1. If there 
exists a standing wave with wave number k in the solid, then 
the displacement of atom m is given by: 



y k (m) — Asm ^/c7r 



n) 



(99) 



The boundary conditions, that ^(0) = 0 and yfe(iV) = 0, 
impose the constraint that k must be an integer, k = 1,2,3.... 
Standing waves with different values of k are known as different 
vibrational modes. 
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Debye model: vibrational modes 
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Debye model: vibrational modes 



From the plot on the previous slide, we see that in mode k = 6, 
the atoms do not vibrate at all: a wave with k = 6 does not 
"exist" as a distinct vibrational mode. 

In fact, further analysis shows that waves with k = 7,8,9... 
simply reproduce vibrational modes with smaller values of k. 
This is because, for example: 

(2N - k)ir— j = -As\n [kTT—J = -y k (m). 

(100) 

The overall phase (minus sign) has no significance for our 
purposes. 
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Debye model: vibrational modes 





Mi 
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Debye model: vibrational modes 



In general, only waves with k = 1, 2, 3, N — 2 exist as distinct 
modes. In other words, there is a maximum wave number, and 
hence a maximum frequency, determined by the number of 
atoms in the solid. 

To determine the maximum frequency, we need to know the 
relationship between the wave number and the frequency. Such 
a relationship is known, in general, as a dispersion relation. For 
a wave propagating through a solid, the dispersion relation 
depends on the forces between the atoms consituting the solid. 

The dispersion relation is the final "ingredient" in the Debye 
model. In the Einstein model, there was only one frequency at 
which the atoms could vibrate. A more realistic assumption, 
often adopted in the Debye model, is that the frequency u of 
the vibration is related to the wave number k by the speed of 
sound c s in the material: 

uj = kc s . (101) 

Statistical Physics 72 Part 5: The Bose-Einstein Distribution 



Debye model: phonons 

Finally, we assume that the energy e in a wave is quantised, 
and that the usual Planck formula applies: 

e = fu J j = hc s k. (102) 

It is important to understand that e is not the vibrational energy 
of a single atom. It represents the energy of a vibrational mode 
in the material, to which many atoms contribute. In essence, 
we are departing from our usual assumption that the energy 
levels in a system of many particles are essentially the same as 
the energy levels of individual particles. To do things properly 
in this case, we should solve the Schrodinger equation for the 
complete system consisting of many particles. However, if there 
are of order 10 23 particles, this is not practical. Instead, we 
arrive at the energy levels using the argument outlined above. 

A vibrational mode with energy e is called a phonon. Phonons 
are quanta of vibrational (sound) energy in a material. 
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Debye model: cut-off frequency 



In three dimensions, the vibrational modes are the same as 
those for photons (electromagnetic radiation) in a closed box. 
The difference is that for phonons, there is a maximum 
frequency (related to the number of particles in the object), 
whereas for photons there is no upper limit on the frequency. 

If g(is)dis is the number of states (vibrational modes) in the 
frequency range v to v + dv, then: 

ru D 



g (v)dv = 3N, (103) 

where vjj is the maximum frequency of a vibrational mode, 
known as the cut-off or Debye frequency. Strictly speaking, the 
right hand side should be less than 3N, since particles at the 
surface of the object cannot vibrate. However, for a 
macroscopic object, the number of particles at the surface is 
small compared to the number of particles in the interior. 
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Debye model: cut-off frequency 



Recall, from Part 3 of this lecture course, that the number of 
states with wave number in the range k to k + dk is: 



g (k) dk = 3x 4tt V ^ k 2 dk, (104) 

(2^)-^ 



where the initial factor 3 accounts for the three "polarisations" 
(two transverse, and one longitudinal) that a phonon can have. 

Using the dispersion relation (101): 

uj = Ikv = kc s , 

we find that the number of states in frequency range v to 
v + dv is: 

V o 

g(v) dv = 12tt^sU Z dv . (105) 

r o 
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Debye model: cut-off frequency 



Combining equations (103) and (105), we find that the cut-off 
frequency is then given by: 

/ g{y)dv=\ 12ir^v 2 dv = 3N. (106) 
Jo Jo ci 



This gives: 



(107) 



The cut-off frequency v D should be the upper limit in any 
integral over the density of states. 
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Debye model: total energy 



Let us proceed to calculate the total vibrational energy U in a 
solid as a function of temperature. Since phonons are bosons 
(they have three possible "spin" or polarisation states), and the 
total number of phonons in a system is not conserved, they 
should follow the same statistical distribution as photons. 

Thus, the total phonon energy is given by: 

U=r D ^de, (108) 

JO e kT — 1 

where e D = hv D is the energy of a phonon at the cut-off 
frequency. 
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Debye model: total energy 



Using equation (104): 

V o 

g (k) dk = 12tt- -=k z dk, 

yy y (2tt) 3 

and equation (102): 

e — hc s k, 

we find that the density of states as a function of energy is 
given by: 

£ 2 

g(e) d£ = 12ttV- -=de. (109) 



Substituting (109) into (108) gives: 



12ttV f£D £ 3 

U = 7T, — \3 L ~ d£ - ^ 11Q ) 



Statistical Physics 78 Part 5: The Bose-Einstein Distribution 



Debye model: total energy 



Then, defining: 

y = —, (HI) 

the expression for the total energy (110) becomes: 

i2irV a r £ D/ kT v 3 
(hcs)^ Jo &y — l 
We define the Debye temperature 0 D : 

s R= hv» 3 = A^ C 3 (n3) 

U k k ' D 4irk 3 V s V J 

In terms of the Debye temperature, the total energy can be 
written: 

T\ 3 i*9 d/t y 3 



u = 9NkT [r D ) Jo ^i dy - (114) 



The definite integral in equation (114) is difficult to evaluate, 
so we consider two limiting cases, of low and high temperature. 
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Debye model: total energy when T<^:6 D 



First, consider the low temperature case, T <.9 D . In this case, 
6 D /T^> 1. Since the integrand becomes extremely small for 
large y, the integral may be approximated by taking the upper 
limit to infinity: 

rO D /T t#3 3 4 

/ — dyn — dy = —. (115) 

JO eV - 1 JO eV - 1 15 

Then, using (114) we find: 

o 4 / rp \ 3 

NkT [ — , \fT<^9 D . (116) 
5 \0dJ 
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Debye model: heat capacity when T <c Ojj 



The heat capacity is given by 

dU 

v 

which gives: 



C v ^ V2 ^Nki [ J^\ , ifT<0 D . (118) 



At low temperature, the heat capacity varies with the cube of 
the temperature: 

Cy 127T 4 Nk 



T 3 5 6% 



as T^O (Debye model). (119) 



This is a different behaviour from that predicted by the 
Einstein model: 

Cy 



T 3 



0 as T^O (Einstein model). (120) 
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The Debye model gives a good prediction of heat capacity 




T 2 \ K 2 ] 



Figure 12.7 — Low temperature molar heat capacity of metals; C Vm /T 
is plotted against T 2 . The results are consistent with equation (12.29), 
according to which C V/n is equal to yT + AT 3 . The intercepts give the 
values of y. The values of A can be determined from the slopes. 
(Reproduced with permission from Statistical Mechanics and Properties 
of Matter by E. S. R. Gopal, Ellis Horwood Ltd., Chichester.) 
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Debye model: total energy when T > flp 



Now, consider the high temperature case, T^>6 D . In this case, 
Qd/T <. 1, and the variable in the integrand remains small 
compared to 1. Then we can approximate e y - 1 « y, and the 
integral becomes: 



o d /t 
0 eV 



V 



dy 



•; !T j*>=\m- 



(121) 



Then, from (114), we find: 

U w 3NkT, 



if T>0 D . 



(122) 
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Debye model: heat capacity when T > 0 D 



Again, we can find the heat capacity: 

At high temperatures, with U = 3NkT (122), we find: 

C v n3Nk, ifT># D . (124) 

At high temperatures, the Debye model and the Einstein model 
give the same value for the heat capacity. 
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Summary: Bose-Einstein Statistics 



You should be able to: 

• Explain that the wave function for a system of identical bosons must be 
symmetric under the interchange of any two bosons; and show that this 
property allows an unlimited number of bosons to occupy any given 
state. 

• Derive the Bose-Einstein distribution for a system of indistinguishable 
bosons. 

• Apply the Bose-Einstein distribution to a Bose-Einstein gas, to show 
the existence of the Bose-Einstein condensation below the Bose 
temperature. 

• Derive the heat capacity of a Bose-Einstein gas, below and above the 
Bose temperature. 

• Explain the difference between a Bose-Einstein gas, and a gas of 
photons. 

• Apply the Bose-Einstein distribution to a gas of photons, to derive 
Planck's radiation law; Wien's law; and Stefan's law (for black body 
radiation). 

• Explain Debye's model of a solid, and discuss the significance of 
vibrational modes and their quantisation (phonons). 

• Use Debye's model to derive expressions for the heat capacity of a solid 
in the limits of low and high temperature. 
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